Universality in ratchets without spatial asymmetry 
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It is demonstrated that to optimally enhance directed transport by symmetry 
breaking of temporal forces there exists a universal force waveform which allows to 
deduce universal scaling laws that explain previous results for a great diversity of 
systems subjected to a standard biharmonic force and provide a universal quantita- 
tive criterion to optimize any application of the ratchet effect induced by symmetry 
breaking of temporal forces. 

PACS numbers: 05.60.-k 

Understanding the ratchet effect [1-4] induced by symmetry breaking of temporal forces 
is a fundamental issue that has remained unresolved for decades. While the dependence of 
the directed transport on each of the ratchet-controlling parameters has been individually 
investigated, there is still no general criterion to apply to the whole set of these parameters 
to optimally control directed transport in general systems without a ratchet potential [5-25]. 
Consider a general deterministic system (classical or quantum, dissipative or non-dissipative, 
one- or multi-dimensional) subjected to a T-periodic zero-mean ac force / (t) where a ratchet 
effect is induced by solely violating temporal symmetries. A popular choice would be the 
simple case of a biharmonic force, fhi,h2(t) = e\har\ {ut + (p x ) + e2har 2 (2ut + (p 2 ), where 
harifl represents indistinctly sin or cos. Clearly, the aforementioned symmetries are solely 
the shift symmetry of the force (/ (t) = —f (t + T/2) , T = 2tt/uj) and the time-reversal sym- 
metry of the system's dynamic equations. Of course, the breaking of the latter symmetry 
implies the breaking of some time- reversal symmetry of the force (/ (— t) = ±/ (t)) in some 
general case, but not in all cases [19]. The analysis of the breaking of these two fundamental 
symmetries allows to find the regions of the parameter space (ei, e 2 , <fi, <f 2 ) > e i + e 2 = const., 
where the ratchet effect is optimal in the sense that the average of relevant observables (such 
as velocity and current, hereafter referred to as (V)) is maximal, the remaining parameters 
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being held constant. In this Letter, it is shown that such regions are those where the ef- 
fective degree of symmetry breaking is maximal. The theory arises from the observation 
that Curie's principle [26] implies that a broken symmetry is a structurally stable situa- 
tion [4]. At this point a quantitative measure of the degree of symmetry breaking (DSB) 
is introduced, on which the strength of directed transport must depend. This quantitative 
relationship between cause (symmetry breaking) and effect (directed transport) is hereafter 
referred to as the DSB mechanism. Also, this quantitative relationship is expected to ex- 
hibit a dependence on the symmetry-breaking parameters which is universal if and only if 
the symmetry breaking takes place solely in the driving force, i.e., in the external agent 
which is simultaneously the transport-inducing force and the ratchet-inducing force. Since 
the ratchet effect can occur at any spatio-temporal scale, such a quantitative measure of 
the DSB must be independent of the force's amplitude and period. I define consistently the 
DSB of the symmetries of the force f(t) by the expressions 

where increasing deviation of -D s ,±(/) from 1 (unbroken symmetry) indicates an increase 
in the DSB. But the effectiveness of any periodic zero-mean force at producing transport 
diminishes as the transmitted impulse over a half-period is decreased while its amplitude and 
period are held constant. In general, this means that optimal enhancement of the ratchet 
effect is achieved when maximal effective symmetry breaking occurs, which is in turn a con- 
sequence of two reshaping-induced competing effects: the increase of the degree of breaking 
of the force's symmetries and the decrease of the transmitted impulse over a half-period, thus 
implying the existence of a universal force waveform which optimally enhances the ratchet 
effect. Thus, for the biharmonic force / cos ,cos (t) = ei cos (uot + y^) + e 2 cos (2ut + </? 2 ), Eq- 
(1) can be put into the form D 8 (f C08jCO8 ) = 1 - (a/n) cos (2r + (p eff ) jP (t; a, (p eff ) dr, 

#+(/cos,cos) = 1 + (°/ 7r ) Io Wsin ( 2r ) Sin( feff/ P ( r ; °> Peff) dl ~ > D ~ (/cos.cos) = 

1 — (1/tt) Jq W [cos r + a cos (2r) cos V e //] /P ( r 5 a > 'Peff) ^ r > where a = e 2 /ei, r = ut+(p 1 , (p = 
(p 2 — 2(p 1 , P (r; a, ^e//) = cos r + a cos (2r + V? e //) ■ ^he quantity (p e ^ is hereafter referred 
to as the effective phase.These integrals diverge at the zeros of the quartic polynomial 
Aa 2 x A + 4acosy?x 3 + (1 — 4a 2 ) x 2 — 2acos(px + a 2 cos 2 y? = 0, where x = cost. After 
solving this algebraic equation for x, one obtains that the three integrals diverge when 
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a ^ 1/2, (p e ff = {7r/2, 37r/2}, and thus the DSB is maximal at these parameter values for 
the three symmetries of the force (see Fig. Id). For these values of the effective phase, 
one finds that the transmitted impulse over a half-period is maximal at a = 1/2 while 
the biharmonic force's amplitude is held constant. This means that maximal effective 
symmetry breaking occurs at a = 1/2, if e ^ = {re/2, 37r/2}. A similar analysis of the 
remaining three versions of the biharmonic force yields the results summarized in Table I 
(second column), which are again the same for the three symmetries in each case. Note 
that one could equivalently define a measure of the DSB by taking the time average of 
the inverse quantities —f(t)/f(t + T/2),±f(t)/f(—t): one finds that the corresponding 
measure (1) exhibits the same qualitative behaviour as a function of the symmetry-breaking 
parameters, and exactly the same optimal values of these parameters are found to yield a 
maximal effective DSB. This indicates that (1) provides a bona fide measure of the DSB. 
Remarkably, such optimal parameter values correspond to a single optimal waveform for 
the four versions of the biharmonic force (see Fig. lb). The DSB mechanism implies that 
such a waveform is universal, i.e., it corresponds to a force waveform which optimally 
enhances the ratchet effect in any system. Consider now the case of the elliptic force 
feiiip(t) — e f{t',T,m,8) = e sn (Qt + 0; m) cn (Qt + O; m), where cn(-;m) and sn(-;m) are 
Jacobian elliptic functions [27] of parameter m, Q = 2K{m)/T, = K(m)9/n, K{m) 
is the complete elliptic integral of the first kind [27], T is the period of the force, and 
9 is the (normalized) initial phase {9 G [0, 27r]). Fixing e, T, and 9, the force waveform 
changes as the shape parameter m varies from to 1 (see Fig. la). In this case, Eq. (1) 
yields D s (f e m p ) = E(m)K~ 1 (m) (1 — m)~ 1 ^ 2 , where E(m) is the complete elliptic integral 
of the second kind [27] (see Fig. lc). Physically, the motivation for this choice is that 
f e iii P (t;T,m = 0,9) = esin (2nt/T + 9) /2, and that f e m p {t;T,m = 1,9) vanishes, i.e., in 
these two limits directed transport is not possible, while it is expected for < m < 1. 
Thus, one may expect in general the average of any relevant observable (V) to exhibit 
an extremum at a certain critical value m = m c as the shape parameter m is varied, the 
remaining parameters being held constant. The DSB mechanism implies that such a value 
m c is universal, i.e., it corresponds to a universal force waveform which optimally enhances 
the ratchet effect in any system. Universality requires that such an optimal waveform 
should be closely related to that deduced for the case of a biharmonic force, in the sense of 
its Fourier series. Indeed, using f e u ip (t;T,m,9)/e = Y^=\ a n(m) sin [n (2nt/T + 9)], 
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a n (m) = mr 2 m 1 K 2 (m) sech [nnK(l — m)/K(m)] one could expect the criti- 
cal value m c to be near m = 0.983417 since f e m p (t;T,m = 0.983417, 9)/e = 
ai (m = 0.983417) [sin (27rf/T + 6) + (1/2) sin (Ant/T + 26) + 0.178592 sin (6vrf/T + 30) + ...], 
i.e., the optimal values (62/61 = 1/2, </? 2 — 2</? 1 =0) for the biharmonic approximation of 
the elliptic function are recovered at m — 0.983417 (cf. Table I, second column, and 
compare Figs, la and lb). Numerical studies of diverse systems [28] confirmed the 
universality and accuracy of the critical value m c = 0.983..., i.e., the universality of 
the optimal waveform. Similarly, from the Fourier series of a sawtooth-wave force 
f sawtooth (t, T) /e = 2 [sin (2id/T) - (1/2) sin (Ant/T) + (1/3) sin (Qnt/T) -...], one recovers 
the optimal values (62/61 = l/2,ty? 2 — 2(/? 1 = ir) for its biharmonic approximation (cf. Table 
I, second column), which explains the great effectiveness of this waveform in controlling 
directed transport of magnetic flux quanta [23]. 

Next, one exploits the aforementioned universality expected from the DSB mechanism 
to deduce the dependence of < V > on the symmetry-breaking parameters (ei,e 2 ,if 1 ,ip 2 ) 
of the biharmonic force fhi,h2(t) in leading order for the usual case [5-22,24,25] of small 
amplitudes (l/ei,2 - ► 00) ■ For the sake of clarity, consider first the case where the viola- 
tion of the time-reversal symmetry of the system's dynamic equations can be absorbed in 
the temporal force because dissipation is negligible and the Lagrangian (Hamiltonian) of 
the system does not contain any additional term explicitly breaking the time-reversal sym- 
metry. This means that the breaking of the time-reversal symmetry implies the breaking 
of the force's symmetry f (—t) = f (t). According to the above arguments, one gener- 
ally expects (V) ~ s (e±, €2, <fi, (p 2 ), where it is assumed without loss of generality that 
the function s is fc-times piecewise continuously differentiable. From MacLaurin's series, 

one has s(e 1 ,e 2 ,<p 1 ,<p 2 ) = Y.T=o Y.n=o C M (<Pu ^2) e i e 2 with c k,o = c , n = since the shift 
symmetry is never broken in the case of a single harmonic function. The transformation 
€i -> -€i,i = 1,2, implies f h i,h2(t) -> -fhi,h2(t), and hence < V >^ - < V >. This means 
that s (ei, 62, (p 1 , <f 2 ) ~ ¥ ~ s ( — e i) — e 2, fi, l P 2 ) and hence k + n = 2m + 1, m — 1, 2, ... . Thus, 
one obtains 

s (ei, e 2 , ifii, <p 2 ) = ci, 2 V^ 2 ) e i e 2 + c 2,i (¥1, <P 2 ) e i e 2 + O {e\e\, e\e 2 2 , e\e 2 , e^el) , (2) 

for ei i2 sufficiently small. Since the ratchet effect does not depend on the time origin, < V > 
must remain invariant under the transformation t — > t + t ,\ft . This transformation yields 
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fhi,hi{t) tihari (cut + (p 1 ) + e 2 har 2 (2ut + </5 2 ), with the fundamental property </5 2 ~~ = 
<f 2 ~ 2ip l , i.e., the effective phase ip 2 — 2<f 1 remains invariant under time translation (see 
Table I, third column), and hence 

(V) ~ ci j2 ((p 2 - 2(pj) t x e\ + C2,i (y? 2 - 2 ^i) efe- (3) 

The transformation ip { — > (p { — n,i = 1,2, implies fhi,hi{t) ~fhi,h2(t), and hence 

< V" >^ - < V > whereby Ci ;2 (</? 2 - 2i Pi) = ~ c i,2 (<f 2 - 2i Pi + , c 2,i (v?2 - 2 ^i) = 
— C2,i (y? 2 — 2(/9 x + 7r), while the transformation ei — > — ei,(/9 1 — > </? 1 — n maintains 
(V) invariant, and hence c\ t2 (<p 2 — 2</? 1 ) = — ci j2 (<f 2 — 2(/? 1 + 2n) , C2,i (</? 2 ~~ 2<^ 1 ) = 
C2,i (<^2 — 2tp 1 + 27r). The comparison of these four relationships for the functions of the ef- 
fective phase implies that c 2j i (<p 2 — 2^) is a 27r-periodic function while Ci j2 (<p 2 — 2<Pi) = 0. 
Thus, Eq. (3) reduces to 

(\/)~(l/ ei )- 2 (l/ e2 )- 1 C2, 1 (^ // ), (4) 

where a power law for the dependence on the amplitudes is now explicit. In this regard, it 
is worth noting the great similarity between the present theory and the highly optimized 
tolerance (HOT) theory [29] where a power law is generated by the actions of an exter- 
nal agent aiming to optimize the behaviour of a system. However, we have seen above 
that universality comes from criticality in the present theory, while for HOT systems the 
details matter. To obtain an explicit expression for the function c (^e//) it * s useful to 
consider the general transformation t — > — t + to,\/t . This transformation yields fhi,H2(t) —* 
e\har\ (out + tp^ + e 2 har 2 (2uut + tp 2 ) where the effective phase is no longer strictly invariant 
but changes according to Table I (fourth column). Since the transformation t — > —t implies 

< V >^ — < V > when the time- reversal symmetry is unbroken, the change rules of the 
effective phase imply that the function c 2j i (f e ff) h &s necessarily a definite parity (cf. Eq. 
(4) and Table I, fourth column). Taking into account this property and given that c 2) i (^e//) 
is assumed to be A;-times piecewise continuously differentiable, its Fourier series [30] can be 
approximated to leading non-trivial order by a single harmonic function according to Table 
I (fifth column). One sees that the universal scaling laws in Table I (fifth column) yield 

< V >= when and only when both the shift symmetry and the time- reversal symmetry of 
the system's dynamic equations (i.e., the force's symmetry / (— t) = f (t) in the present case) 
are unbroken, while they yield a maximum value of < V > when and only when maximal 
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effective symmetry breaking occurs in the sense of the measure (1), as predicted from the 
DSB mechanism. Also, that the harmonic functions appearing in the universal scaling laws 
are independent of har\ is a consequence of the invariance of < V > under time translation. 
As expected, one finds that such universal scaling laws confirm and explain previous results 
for a great diversity of systems [7,14,15,22,25] subjected to a biharmonic force fhi,h2(t)- 

I now discuss how the aforementioned universal scaling laws change when the violation 
of the time-reversal symmetry of the system's dynamic equations cannot be absorbed in 
the temporal force. This is the case when dissipation [5,6,8-13,16-21,24] is not negligible. 
It has been demonstrated above the approximate conservation of the effective phase in the 
sense of the change rules in Table I (fourth column), and hence that ip e ^ is the proper ar- 
gument of the function C2,i (tp 1 , tp 2 ), when the violation of the time-reversal symmetry only 
occurs in the temporal force. Therefore, that the violation of such a symmetry is also due 
to the presence of dissipation means that ip e ff can no longer be an argument of the function 
C2,i (<£>i, (p 2 ) but one has ip e ff + (p diss instead, where (fi diss is hereafter referred to as the dis- 
sipation phase. Note that the additive character of the dissipation phase is a consequence 
of the DSB mechanism. Thus, the dissipation phase quantifies the degree of breaking of 
the time-reversal symmetry generated by dissipation. Also, the DSB mechanism implies 
the universal properties: (p diss (f3 — 0) = and that (p diss (f3) is a monotonously increasing 
function of (3, with (3 being the effective dissipation parameter. For the values of y? e // yield- 
ing < V >= in the absence of dissipation, i.e., those values for which the temporal force 
does not break the time-reversal symmetry of a non-dissipative system, one obtains that 
the maximum absolute value (i.e., 1) of the harmonic functions appearing in the universal 
scaling laws is reached at tp diss = ±vr/2, and hence we have the additional universal property 
ma.x /3 ip diss (P) = n/2 (cf. Table I, fifth column). However, the function (p diss (/3) generally 
depends upon additional parameters, such as the period and diverse system-dependent pa- 
rameters, i.e., it is not a universal function. Of course, it is generally expected that the 
function < V > jhar (v e // + Pdiss) should exhibit monotonously decreasing behaviour as a 
function of /3, where har is the corresponding harmonic function in Table I (fifth column) 
in each case. When dissipative forces dominate inertia (the so-called overdamped regime 
[1,4]), the breaking of the time- reversal symmetry implies the breaking of the force's symme- 
try / (— t) = — fit) and the dissipative phase reaches its limiting values f diss = ±7r/2. Since 
the optimal values of the relative amplitude e 2 /ei and the effective phase ip e fj are just the 
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same for the three symmetries of the biharmonic force, this means that the universal scaling 
laws corresponding to the overdamped regime are those given in Table I (fifth column) but 
with sin instead of cos, and vice versa, in each case. One finds that these predictions are in 
perfect agreement with published results for a great diversity of systems [4-6,8-13,16-21,24] 
subjected to a biharmonic force fhi,h2(t). Since dissipative forces and randomly fluctuating 
forces (noise) have the same microscopic origin, it is expected the effectiveness of temporal 
forces at generating directed transport induced by the ratchet effect to be robust against 
moderate presence of noise. 

In summary, universal scaling laws for the strength of directed transport induced by 
symmetry breaking of temporal forces have been deduced from a quantitative interpretation 
of Curie's principle. The present theory explains in a general setting all previously published 
results for a great diversity of systems [5-25], and provides a universal quantitative criterion 
to optimize any application of the ratchet effect induced by symmetry breaking of temporal 
forces. 
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TABLE I 



hari, har 2 


D s ,± 


t -> t + t 




(V) 


cos, COS 


U = V2, ¥>«// = {§,¥} 


= ¥>e// 


Veff = ~<Peff 


~ e\e 2 sm(p eff 


sin, sin 


a = 1/2, v, e// = {0, 7T> 




<Peff = -<Peff ± * 


~ e\e 2 cos y? e// 


sin, cos 


S = l/2^ e// = {f,f} 


= ¥>e// 


<Peff = -<Peff 




cos, sin 


a = l/2^ e// = {0,7r} 


$eff = <Peff 


<Peff = ~<Peff ± * 


~ efe 2 cos y? e// 



Table I. Optimal values of the relative amplitude £2/61 and the effective phase ip e ^ = 
ip 2 — 2(/? 1 obtained by computing the measure (1) of DSB for the three symmetries of the 
biharmonic force (second column), change rules of the effective phase under time transfor- 
mations (third and fourth columns), and universal scaling laws in leading order for averaged 
velocities and currents. Note the coherence of the results in the second and fifth columns, 
which were obtained using independent methods (recall that, without loss of generality, 
ei + e 2 = 1 so that e\e 2 = (1 — e 2 ) 2 e 2 , which is a function having a single maximum at 
e 2 = 1/3, and hence e\ = 2/3,e 2 /e 1 = 1/2.) 

Figure Captions 

Figure 1. (a) Elliptic force f e iu P (t) = ef(t; T, m,9) = e sn (Qt + 0; m) cn (fit + 0; m) vs 
t/T and three shape parameter values, m = (light blue line), m = 0.983417 (dark blue line, 
optimal universal waveform), m = 1 — 10 -6 (light blue line), showing an increasing symmetry- 
breaking sequence as the pulse narrows, i.e., as m — > 1. (b) Optimal universal biharmonic 
force generating directed transport in one direction (f^ har (t) je — sin(27rf/T) + | sin (Ant/T), 
dark blue line) and the opposite direction (f^ har (t) /e = sin(27rt/T) — (1/2) sin(47rt/T) = 
—fbih ar {t-\-T/2) /e, red line), (c) Measure of the DSB (Eq. (1)) for the elliptic force in 
(a), D s (f e m p ) = E(m)K~ l (m) (1 — m) -1 ^ 2 vs m. One sees a sharp increase as m — > 1. 
(d) Measure of the DSB (Eq.(l)) for the biharmonic force fbihar(t)/e = cos(27rt/T) + 
a cos (4-7rt/T + <f e ff), D s (f bihar ) vs a for ip e ^ = ir/2. One sees a sharp increase as a — > 1/2, 
which is similar to that found in (c) for the elliptic function. 



